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ABSTRACT 

Abundances of N, Si, S, and Fe for 45 damped Lyman alpha systems (DLAs) 
have been compiled and detailed one-zone chemical evolution models have been 
constructed for 30 of them. Assuming continuous star formation, we found that 
final abundances in each object can be modelled by adjusting only two param- 
eters, i.e. its time-averaged star formation efficiency and evolutionary age, with 
ranges in our sample of 0.01-1.5 Gyr -1 and 0.18-2.0 Gyr, respectively. In addi- 
tion, average star formation efficiency and evolutionary age appear to be anticor- 
related for the sample, suggesting that the star formation efficiency in a typical 
DLA decreases with age. At the same time, N/Si in DLAs is directly linked 
to an object's age. There is an apparent bimodality in the distribution of N/Si 
values which could be the result of a statistical accident or an effect produced 
by a truncated or flattened IMF. We find that the mean and small dispersion 
of Si/Fe values is related to the generally young ages of DLAs, wherein not all 
Fe has yet been expelled by Type la supernovae. Finally, the large scatter and 
generally lower values of N/Si of DLAs with respect to blue compact galaxies, 
despite their partially overlapping metallicities, indicate that DLAs are generally 
younger than the latter. 
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1. INTRODUCTION 

Damped Lya systems (DLAs) serve as valuable probes of heavy element abundances and 
chemical evolution at high redshifts and thus large lookback times (Lu et al. 1996; Pettini, 
Lipman, & Hunstead 1995). They provide complementary information to that obtained from 
stellar and nebular studies of nearby galaxies (Pritzl, Venn, & Irwin 2005; Bresolin, Garnett, 
& Kennicutt 2004). 

The canonical view of DLAs is that they are galactic systems, perhaps representing a 
variety of morphological types, whose interstellar matter absorbs light from QSOs positioned 
behind them along the observer's line of sight. These systems have H I surface densities 
where N(H I) > 2.0 x 10 20 cm~ 2 (roughly 2 M©pc~ 2 ), i.e. high enough so that absorption 
occurs in both the core as well as the damping wings of the Lja feature. Under such 
optical depth conditions metal lines become detectable, and careful analysis of their profiles 
enables researchers to measure abundances of more than 20 elements in some DLAs (e.g. 
Prochaska, Howk, & Wolfe 2003). For a comprehensive discussion of numerous aspects of 
DLAs, including their contributions to the cosmic abundance picture, the reader is referred 
to the recent review by Wolfe, Gawiser, & Prochaska (2005). 

Two approaches to understanding DLA abundances exist in the literature. The "global" 
one treats abundance averages over co-moving volumes and attempts to model the observa- 
tions using global values of star formation histories (Edmunds & Phillipps 1997). In contrast, 
the "local" approach studies individual DLAs one at a time, using standard chemical evolu- 
tion equations to interpret observations for a sample of objects (Dessauges-Zavadsky et al. 
2006, 2007). We focus on the local method here. 

Confronting abundance measurements in DLAs with numerical models has proved fruit- 
ful for understanding the chemical evolution of these objects. Using this approach, Matteucci, 
Molaro, & Vladilo (1997) found that abundance patterns in DLAs are consistent with burst- 
ing star formation in a dwarf irregular galaxy, thus implicating the latter as progenitors for 
some DLAs. Later, a similar analysis by Calura, Matteucci, & Vladilo (2003) appeared to 
rule out elliptical galaxies as possible progenitors but could not decide between spiral disks, 
irregular galaxies, or bursting dwarf irregulars for the host. Dessauges-Zavadsky et al. (2004) 
employed both a spiral disk and dwarf irregular galaxy model in efforts to explain abundance 
measurements in three objects and concluded that DLAs are either outer regions of spiral 
disks which have undergone slow continuous star formation or are dwarf irregular galaxies 
which have experienced either continuous or bursting star formation. Their models had the 
added virtue of predicting star formation rates which were consistent with the rates inferred 
by Wolfe, Prochaska, & Gawiser (2003) from observed C II* absorption. The recent paper 
by Dessauges-Zavadsky et al. (2006) used models from Dessauges-Zavadsky et al. (2004) 
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but extended their survey to report abundances of 22 elements for a sample of 11 objects. 
Finally, Dessauges-Zavadsky et al. (2007) has computed chemical evolution models of six 
additional DLAs using the same approach as they did in the 2004 paper and reached similar 
conclusions regarding the plausible star formation histories for these objects. 

Metallicity in DLAs [M/H], as traced by [Si/H], [S/H], and [Zn/H0, ranges from roughly 
—0.5 to — 2.6dex , with a cosmic mean metallicity [M/H]« — 1.1 dex (Prochaska & Wolfe 
2000; Prochaska et al. 2003; Pettini 2004). The observations indicate the metallicity increases 
with the age of the universe and that no DLA has a metallicity extending below solar by 
more than -2.6 dex (Prochaska et al. 2003). These conclusions have been bourne out by more 
recent compilations (Meiring et al. 2006; Ledoux et al. 2006; Prochaska et al. 2007). The 
typical values of the alpha abundance, as traced by Si and S, are 1-2 orders of magnitude 
below solar values. 

DLA abundance studies over the past decade have focused much of their attention on 
two important element /element ratios, i.e. N/a and a/Fe, where the relevant alpha element 
is usually Si, with O and S sometimes playing the role. The choice of Si is not ideal given 
that this element is mildly refractory, but the resonance- line O I transitions are generally 
saturated and the triplet of S II transitions near A « 1250A are frequently lost in the Lya 
forest. 

The characteristic range of the N/a ratio, as well as the manner in which objects appear 
to form a bimodal distribution, has sparked much debate. These issues are of particular 
interest in determining overall N production in galaxies, because the evolution of this ratio is 
linked to the relative roles of both low and intermediate mass stars (LIMS) on the one hand 
and massive stars on the other. Pettini et al. (2002) employed VLT UVES observations 
along with published data to study the N/O ratio in 10 DLAs and confirmed the broad 
range in values seen earlier by Lu et al. (1998). Pettini et al. proposed that this dispersion 
could be explained if LIMS produced a large fraction of the N, and that their sample objects 
represented various stages of N release by LIMS, which evolve somewhat slower than massive 
stars. This proposal was consistent with chemical evolution models of Henry, Edmunds, & 
Koppen (2000; hereafter HEK). Then Prochaska et al. (2002) reported a possible bimodality 
in the distribution of N/a later supported by Centurion et al. (2003). They found that while 
most DLAs show [N/a] m —0.5, a small group of objects (dubbed "low nitrogen DLAs" 
or LNDLAs) have [N/a] ~ —1.5. They considered several explanations for this bimodality, 
including the delay in N release favored by Pettini et al. and concluded that this would not 
likely result in a bimodality, claiming that under this scenario the development of a potential 



1 We adopt the standard bracket notation, [X/H] = log(X/H) — log(X/H)0 
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gap between the normal N objects and the LNDLAs would likely not occur. Instead, their 
preferred hypothesis was a reduction in the effective N yield effected by a highly flattened 
or truncated IMF possessing a deficit of LIMS. (Prochaska et al. assumed that LIMS were 
responsible for most of the N production in the Universe, an issue which is yet to be settled.) 

Following the claim of an N/a distribution bimodality, Molaro (2003) and D'Odorico 
& Molaro (2004) confirmed the bimodality and argued in favor of the delay hypothesis. 
Lanfranchi & Matteucci (2003) showed with detailed models that the top-heavy IMF proposal 
of Prochaska et al. (2002) could indeed explain the low N/a values in LNDLAs, but at 
the expense of significantly underproducing Fe, with the result that the predicted Si/Fe 
ratio is much larger than observed. More recently, Chiappini, Matteucci, & Meynet (2003) 
and Chiappini, Matteucci, & Ballero (2005) have employed the stellar yields of Meynet & 
Maeder (2002) to show that the LNDLAs could arise in systems undergoing short bursts (10- 
200 Myr), while the star formation history of most DLAs is characterized by significantly 
longer bursts (800 Myr). Finally, the N/a problem was addressed in the work of Molla et 
al. (2006), who showed that by using the yields for LIMS computed by Buell (1997) and 
subsequently published and used in Gavilan, Buell, & Molla (2005) and Gavilan, Molla, & 
Buell (2006), they could explain a large range in N/a in astronomical objects, including the 
high values recently found in halo stars by Spite et al. (2005) and Israelian et al. (2004), the 
moderate values seen in blue compact galaxies, and normal as well as low nitrogen DLAs. 
Clearly, the problem of N/a is far from resolved. 

Another major problem concerning DLAs is the attempt to explain the supersolar a/Fe 
ratios observed in nearly all DLAs. For example, in their extensive study of DLA abundances 
Prochaska & Wolfe (2002) determined a mean value of log(Si/Fe) of +0.3 dex with a small 
standard deviation. This finding is consistent with the recent results of Dessauges-Zavadsky 
et al. (2006) who derived a mean value of log(Si/Fe) of +0.430 dex with a standard deviation 
of 0.114 dex. The majority of authors (e.g. Vladilo 1999; Pettini et al. 2000; Centurion et al. 
2000; Nissen et al. 2004) have argued that the observed supersolar ratios of Si/Fe are simply 
the result of depletion; if the Fe abundance is corrected for depletion the Si/Fe becomes solar. 
However, Wolfe, Gawiser, & Prochaska (2005) plot [Si/Fe] versus [Si/H] and argue that in 
the limit of vanishing metallicity, where dust depletion is insignificant, [Si/Fe] approaches 
~ +0.3, not zero. Furthermore, Dessauges-Zavadsky et al. (2006) have fitted the Si/Fe 
ratio as a function of observed Zn/Fe value (a measure of differential dust depletion) and 
found supersolar Si/Fe ratios for "dust-free" systems. Resolution of the Si/Fe ambiguity 
is of crucial importance for understanding the chemical nature of DLAs; a true solar value 
implies that these systems are much more advanced chemically than if Si/Fe is significantly 
above solar. 
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The goal of this paper is to explore the N/a and a/Fe problems just discussed for DLAs, 
as well as other questions related to their metallicity range, by employing a very basic chem- 
ical evolution model with few a priori assumptions along with the smallest number of variable 
parameters. To this end we have compiled published abundances of N, Si, S, and Fe for a 
sample of 45 DLAs. Using only those objects for which firm values of N and Si abundances 
are available (that is, objects with only limits for one or both of these elements were re- 
jected), we computed a detailed chemical evolution model for each object, where the only 
parameters which were varied were the star formation efficiencj|§ (SFE) and evolutionary 
age (measured from the time when star formation began to the observed epoch.). We then 
used the results to address the abundance problems presented above. 

Section 2 presents and discusses the abundance compilation upon which our analysis is 
based, while § 3 describes the implications of the data. The results and implications of the 
chemical evolution models are provided in § 4, and our conclusions are presented in § 5. 

2. THE DATABASE 

Table CD provides a summary of all of the DLAs with N constraints in the literature. 
Columns 1 and 2 contain the name of the background quasar and redshift of the absorbing 
system, while columns 3-10 give abundance information. The literature sources for the 
abundances are given in the last column. 

All DLAs were observed with either the HIRES (Vogt et al. 1994) or UVES (Dekker 
et al. 2000) echelle spectrometers on the Keck 10m and VLT 8m telescopes, respectively. 
The data have high resolution (FWHM < lOkms -1 ) and moderate to high signal-to-noise 
(S/N > 10pix _1 ). The sample includes ~ 10 new observations acquired with the HIRES 
spectrometer using the upgraded mosaic CCD (Prochaska et al. 2007, in prep). These data 
have been reduced and analyzed with the same techniques described in Prochaska et al. 
(2002). 

The nitrogen and a elemental abundances were derived from low-ion transitions and 
reflect gas-phase abundances. The ionic column densities were measured with line-profile 
fits or the apparent optical depth method (Savage & Sembach 1991). The reported errors 
only indicate statistical uncertainty. One may allow for an additional ~ 0.05 dex error 
for systematic uncertainty, e.g. continuum placement. This is especially true for the NI 



2 The star formation efficiency (SFE) is the mass of new stars formed per mass of material available for 
forming stars per unit time. The unit for the SFE is therefore Gyr -1 . 
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transitions, all of which lie within the Lya forest. 

In general, one must consider ionization corrections and differential depletion in order to 
convert the gas-phase ionic column densities to intrinsic elemental abundances. Conveniently, 
N, S, and O are non-refractory and Si appears to be only mildly depleted (Savage & Sembach 
1996; Welty et al. 2006); [Si/S] ~ in those DLA where both are measured (Prochaska & 
Wolfe 2002; Dessauges-Zavadsky et al. 2006). Ionization corrections, however, may be 
important for nitrogen, which has a relatively large cross-section to X-rays (Sofia & Jenkins 
1998). Approximately half of the DLAs listed in Table [TJ have had their ionization states 
assessed in detail, and several DLAs have been corrected by a few 0.1 dex, while others have 
been removed from further analysis because of large uncertainty (Prochaska et al. 2002). We 
suspect that a few of the other DLAs drawn from the literature, e.g. Centurion et al. (2003), 
may require modest ionization corrections which would increase the N abundance by 0.1 to 
0.3 dex. While these systematic corrections exceed the reported statistical uncertainty, we 
emphasize that effects of this order have little bearing on our conclusions. 

Observed values of log(N/Si) and log(Si/H) from Tabled] are plotted against each other 
in Fig. [TJ The filled circles show 30 points for DLAs with firm values, while up/down open 
triangles indicate an additional 15 objects where only lower/upper limits for log(N/Si) could 
be determined. For comparison purposes we also show positions for 68 blue compact galaxies 
(BCG) taken from Nava et al.'s (2006) nitrogen study. In order for the BCGs to be placed 
on this plot the O abundances given by Nava et al. were converted to Si abundances using 
log(Si/O) = -1.15 (Asplund et al. 2005). We remind the reader that each object represented 
in Fig. [1] has evolved independently from the others and has a unique star formation history. 
Each point indicates an object's current position in N/Si-Si/H space. 

It is interesting to note the relatively large range in log(N/Si) exhibited by DLAs com- 
pared with BCGs. Nava et al. defined an N/O plateau as comprising those objects forming 
a narrow horizontal track whose vertical dimension is indicated with a double-headed arrow 
in Fig. [TJ We note that all but a few DLAs with firm values for N/Si lie on or below this 
plateau, with Q1008+36 appearing to have the smallest N/Si ratio, -1.15, of all the DLAs 
in Table [TJ In contrast, the DLA with the highest value of log(N/Si), +0.22, is J1558-0031, 
a system with relatively simple kinematics that also shows D absorption (O'Meara et al. 
2006). This vertical range of -1.15 to +0.22 for DLAs can be compared with the range in 
N/O for the plateau of BCGs given by Nava et al. (and converted to N/Si) of -0.39 to -0.12. 
The log of the mean N/Si values for DLAs is -0.40 (-0.46 when J1558-0031 is excluded), 
compared to -0.28 for BCGs. Restricting the BCG sample to the —0.39 to —0.12 interval 
and the DLA sample to log(Si/H) > —6, a two-sided Kolmogorov-Smirnov test rules out 
the null hypothesis that the N/Si distributions of the two samples are drawn from the same 
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parent population. Finally, Fig. [T] shows that DLA metallicities, as measured here by Si, 
extend below the lowest BCG metallicity by as much as 0.6 dex. 

The LNDLAs discussed in the Introduction are visible as the group of four objects with 
firm abundances (filled circles) and one upper limit (open triangle) located in the lower left 
of Fig. [1] inside the box. Focusing on the circles for a moment, there is an apparent void 
stretching obliquely from upper left to lower right that separates these objects from the bulk 
of the DLA population. 

A second interesting pair of quantities to consider is log(Si/Fe) versus log(Si/H), as the 
first parameter illustrates the a/Ye characteristics of DLAs. Only objects in Tabled] with 
firmly measured abundances of Si and Fe are shown in Fig. [2l The solar value for log(Si/Fe) 
from Asplund et al (2005) is +0.06 and is indicated with a horizontal dashed line. The 
vertical range in log(Si/Fe) extends from +0.13 to +0.89, while the log of the mean of Si/Fe 
is +0.44. With a solar value for log(Si/H) -4.49, the DLAs fall roughly between -0.71 and 
-2.3 dex below solar metallicity. Our figure should be compared with a similar one in the 
bottom panel of Fig. 8 in Wolfe, Gawiser, & Prochaska (2005) in which it is suggested that 
Si/Fe is constant at low Si/H but turns upward when [Si/H]>-1. The general consensus is 
that the increase in Si/Fe above one-tenth solar metallicity is due to increasing depletion of 
Fe onto dust as the metallicity rises (Prochaska & Wolfe 2002). 

This general pattern is very similar to what is observed for stars of low metallicity in 
the Milky Way Galaxy and dwarf spheroidals. Recently, Pritzl, Venn, & Irwin (2005) have 
compiled numerous abundances for Galactic globular cluster stars. Data for their objects 
which fall within the Si/H range of our plot are shown as crosses in Fig. [2j Keeping in 
mind the likely differences in the star formation histories of DLAs and the Milky Way, we 
nevertheless find that the log of their mean value for Si/Fe for all of their objects is +0.44, 
identical to our DLA value. Supersolar values of Si/Fe are thought to arise in stars when 
they form from material which has been enriched primarily by massive star production of 
alpha elements such as Si, while type la supernovae have yet to expel completely their Fe 
peak products. 

Finally, in Fig. Owe present a plot of log(S/H) versus log(Si/H) for the 19 objects in 
Table [1] with firm values for both of these elements. The bold line shows a least squares fit 
to the data, and the elements of the fit are given in the upper left corner of the graph. The 
correlation coefficient (r) and slope (B) are both near unity, signifying a solid linear relation 
between these two elements. The value for A, -0.31, is nearly equal to the solar value for 
log(S/Si) of -0.37 given by Asplund et al. (2005). 

In summary, the following points can be made about the patterns seen in the data: 
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1. The N/Si ratios exhibit a large amount of scatter for DLAs compared with the relatively 
narrow distribution in N/O observed in blue compact galaxies. 

2. N/Si and Si/H in DLAs generally fall below corresponding values for blue compact 
galaxies. 

3. The distribution of DLAs in the N/Si-Si/H plane appears to be bimodal (see § 4.2.1). 

4. Observed values of Si/Fe for DLAs are consistently above the solar level for the same 
ratio. 

5. A tight correlation exists between S/H and Si/H. 

We now endeavor to explore the above characteristics of DLAs using chemical evolution 
models. In doing so we focus only on those objects for which firm values for both N and Si 
abundances are available; objects for which only limits have been determined for either N or 
Si are ignored. This culling reduced our sample to 30 objects. 

3. ANALYSIS 

An individual chemical evolution model for each of the 30 sample objects was calculated, 
after which model predictions of N/Si, Si/Fe, Si/H, S/H, and Fe/H were compared with 
their observed counterparts. Details of the chemical evolution code that we employed are 
provided in Appendix A. Each model was a one-zone calculation, that assumed continuous 
star formation, the presence of infall (but not outflow), a stellar mass range between 0.1 and 
120 M , a Salpeter (1955) IMF, no instantaneous recycling, instantaneous mixing, and a 
standard set of stellar yields. This is our basic model. The only quantities which varied from 
model to model were the input SFE and evolutionary age. The abundance levels of each of 
the isotopes 4 He, 12 C, 14 N, 16 0, 20 Ne, 28 Si, 32 S, 35 C1, 36 Ar, and 56 Fe were tracked over time. 

For our entire analysis, massive star yields were taken from Portinari et al. (1998) for all 
of the above isotopes but 35 CI and 36 Ar; yields for these last two species came from Woosley 
& Weaver (1995). Yields of LIMS for 4 He, 12 C, 14 N came from Marigo (2001). Neither yield 
set accounts for stellar rotation. Contributions of SNIa to the evolution of all but 4 He were 
included using the yields for the W7 and W70 models of Nomoto et al. (1997). All carbon 
yields were multiplied by 0.3 in keeping with our findings in Henry, Nava, & Prochaska 
(2006). 

Table [2] shows the results of the model of each DLA in our culled sample of 30 objects. 
Object names are provided in column 1, while columns 2 and 3 list the SFE and age used 
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as input in the model. Columns 4-12 list the observed and model predicted values for the 
quantity specified in the column heading. To facilitate comparisons, each table element is 
given in the format observed/model. Ellipses are entered when no observational values are 
available. Since few DLA carbon abundances have been determined, only model predictions 
are given in columns 9 and 12. The last two rows of the table contain the average observed 
uncertainty relevant to each column and the solar abundances from Asplund et al. (2005), 
respectively. 

During the modeling process, the principal goal was to adjust the SFE and age in order 
to reproduce specifically the observed values of N/Si and Si/H. Therefore, observation and 
theory are in particularly close agreement for these two ratios (columns 10 and 5, respec- 
tively). Sulfur and iron predictions, while not primary modeling targets, are nevertheless 
in good agreement with the observations. Major discrepancies occur in the case of oxygen, 
however, where typically the models overpredict the abundance by an average of 0.43 dex. 
We speculate that part of this difference may be due to line saturation in the observed O I 
transitions, causing observed O/H to be underestimated. 

Only a small number of our objects have C and O abundances available in the literature. 
DLA carbon abundances in particular are extremely difficult to measure because the only 
resonance line outside of the Lya forest, C II A1334.5, becomes saturated even in metal-poor 
systems. We find that the ratio of C/Si is roughly the same throughout our DLA sample, 
with an average value in our models of +0.22 (er=.059). D'Odorico & Molaro (2004) observed 
a value of +0.66 dex for C/Si in BR1202-07. Centurion et al. (2003) report O abundances for 
three of the objects in our sample, and these observations are included in Table [2j We predict 
an O abundance only 0.04 dex larger than the observed value for Q0000-2619, although our 
predictions for Q0347-38 and Q2206-19 are 0.23 dex and 0.33 dex less than measured levels, 
respectively. We also found that O/H scales tightly with Si/H in our models, as expected 
from massive star nucleosynthesis. 

Fig. H] shows plots of residuals between model predictions and observed abundances 
versus the corresponding observed values (i.e. no dust corrections) for N/H, Si/H, Fe/H, and 
S/H. The mean of the absolute values of the residuals is also indicated in each panel. Dashed 
lines indicate levels of zero residual. Again, we see that agreement between observations and 
the model predictions is very good in particular for N/H, Si/H, and S/H, where the average 
residual is less than 0.1 dex in each case. The average residual for Fe/H is only slightly 
larger. We also see a preference for model predictions of S/H to exceed observations, with 
an even larger trend in this direction in the case of Fe, which may imply mild depletion of 
Fe onto dust. 

Finally, the uniqueness of models in Table [2] was assessed using Monte Carlo simulation 
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techniques in conjunction with the parameter ranges inferred from our models. The test 
is described in detail in Appendix B, where we conclude that, within the confines of our 
assumptions, our models have a high probability of being unique. 

4. DISCUSSION 

4.1. Star Formation Efficiency and Evolutionary Age 

The values for log(SFE) in Gyr -1 and log(age) in Gyr for each object, as inferred from 
our models, are plotted against each other in Fig. [51 The four LNDLAs are enclosed in 
a box. Note that the LNDLAs are grouped together in the region of low SFE and age, 
consistent with their positions in Fig. [T] at low values of Ni/Si and Si/H. Values for SFE 
range between 0.01-1.5 Gyr -1 , with a mean and median of 0.44 and 0.20 Gyr -1 , respectively. 
At the same time, the evolutionary ages range between 0.18-2.2 Gyr. The SFE values for 
DLAs are similar to those inferred for BCGs (Henry, Nava, & Prochaska 2006), although 
BCGs tend to be more metal-rich (see our discussion of Fig. [TJ above). They are also in 
good agreement with the model results of Lanfranchi & Matteucci (2003) and Chiappini, 
Matteucci, & Meynet (2003) for BCGs and DLAs, suggesting that BCGs are more advanced 
than DLAs in the evolutionary sense. Clearly the DLAs in our sample cannot have evolved 
for a time greater than the ages implied by the redshifts of the absorbing systems. At the 
same time, evidence is building that indicates that BCGs contain AGB stars and are at least 
several Gyr in age. (See Aloisi, Tosi, & Greggio 1999 and Ostlin 2000 for the specific case 
of I Zw 18). 

One of the most striking attributes of Fig. [5] is the apparent anticorrelation between SFE 
and age, which is especially strong among the objects not enclosed in the box, i.e. those 
which are not LNDLAs. The SFE inferred through our models is actually the time- averaged 
value over the age of a system, since the SFE remains constant throughout the calculation. 
But star formation could actually have been much more prodigious early on in many objects. 
This would be the case if bursts dominated at an early age, and then gave way to a more 
tranquil period of quiescent, continuous star formation. Supposing then that all DLAs in 
our sample evolved in this way, a younger object would have a higher average SFE than 
an older one. Ultimately, this finding could be related to the decline in surface density of 
star-forming materials with time as a system forms stars and evolves. We conclude that the 
SFE declines as a DLA evolves. 

At the same time, the fact that we derive a supersolar Si/Fe value for all of our objects 
suggests that star formation has not ceased altogether. If it had, Si production by rapidly 
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evolving massive stars would have stopped while Fe production from slower evolving Type la 
supernova progenitor systems would have continued, resulting in more solar-like values for 
Si/Fe (e.g. Dessauges-Zavadsky et al. 2004). 

Finally, the data in Fig.[5]imply that there are no objects characterized by both relatively 
high SFE and advanced age, suggesting that DLAs with those characteristics are either very 
rare or difficult to detect. 

In relatively unevolved objects such as DLAs where the gas fraction is nearly unity 
(see eq. A10), the star formation rate is equal to the product of the SFE and the surface 
mass density, i.e. ip ~ uM. While our models help us constrain the SFE, they do not limit 
the surface mass density, since our predicted element-to-element ratios are not a sensitive 
function of that quantity. Thus, we are unable to use our results to infer unique star formation 
rates to compare with observed ones. However, by combining our derived SFE values with 
the average DLA star formation rates from Wolfe, Prochaska, & Gawiser (2003) we can infer 
something about the surface densities of the DLAs within our sample. 

For an average star formation rate Wolfe et al. found log^=-2.2 when a cold neutral 
medium (CNM) is assumed, and log^—1.3 for a warm dark medium (WNM), where the units 
of ip are M Q /kpc 2 /yr. In the case of a CNM then, we find a range of surface mass densities 
in M Q /pc 2 of 3.2 to 315.5 with average and median values of 69.1 and 31.6, respectively. 
Likewise, for the WNM scenario, we find the range to be 25.05-2505 with an average and 
median of 548.6 and 250.5. While the WNM results appear to give unrealistic values for 
surface mass density, the CNM results are for the most part consistent with typical values 
of dwarf irregulars (the low end of the density range) and spiral disks. 

Correlations between abundance ratios on the one hand and SFE and evolutionary age 
on the other are important for understanding the underlying nature of DLAs. In Fig. [6] we 
plot observed values of log(N/Si) and log(Si/H), in the upper and lower panels, respectively, 
versus model- determined values of log(SFE), log(age), and log(SFE x age) left to right, 
respectively, as independent variables. Note that log(SFE x age) is related to total astration. 
In the upper panels, there appears to be a weak anticorrelation between log(N/Si) and 
log(SFE), a very strong direct one between log(N/Si) and age, and no apparent relation 
between log(N/Si) and log(SFE x age). Note also that the strong link between N/Si and 
age that we find is consistent with the suggestion by Edmunds & Pagel (1978) regarding 
extragalactic H II regions that N/O tracks age and is due to the sensitivity of N evolution to 
the delayed release of this element by LIMS as they evolve at a slower pace relative to massive 
stars. The weak anticorrelation between N/Si and SFE is linked to the anticorrelation 
described above between SFE and age. The lower three panels show a weak direct link 
between Si/H and SFE but a very strong one between Si/H and SFE x age. Since the latter 
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quantity is related to total astration (the total number of stars formed over the lifetime of 
the system) and manifested in a system's metallicity, it is not surprising that SFExage is 
closely linked to a metallicity tracer such as Si/H. 

In computing our models, then, we relied upon the tight correlation between log(N/Si) 
and age to constrain the latter parameter. Once the age was determined, we could then adjust 
the SFE to bring the predicted value of log(Si/H) into line with its observed counterpart. 
Thus, despite the apparent degeneracy between SFE and age demonstrated in Fig. we 
were able to infer a unique set of input parameters for each object in this way. 

Only a small number of chemical evolution models of DLAs have been published, so a 
comparison with others' results is very limited. Dessauges-Zavadsky et al. (2004) modeled 
three objects, including Q1331-17, which is also present in our sample. These authors cal- 
culated both a spiral galaxy and a dwarf galaxy model, with the latter being more directly 
comparable to ours for this object. They derived a value for SFE of 0.03 Gyr -1 , compared 
with our finding of 0.06, along with an age > 3.5 Gyr compared to our value of 1.2 Gyr. 
Their greater age was necessary to match the near-solar Si/Fe value which they inferred 
after applying a substantial dust correction to the observed Si/Fe ratio. 

In a more recent paper which appeared after our study was nearly complete, Dessauges- 
Zavadsky et al. (2007) used the same approach as in their 2004 paper to model six more ob- 
jects, four of which are present in our list, i.e. Q0841+129 (z=2.375), Q1210+17, Q2230+02, 
and Q2348-1444. Comparing our SFE and age results with the same parameters from their 
favored models the two studies are consistent for all but the case of Q2230+02, where we 
derive a much larger SFE and much smaller age. 

Finally, whatever the details are regarding the physical systems hosting DLAs, their 
chemical evolution is currently comprehensible simply through the use of very basic one- 
zone models characterized by low continuous star formation rates and only two variable 
parameters, namely the system's SFE and evolutionary age. At this point in our under- 
standing of these objects, no other specifications appear to be needed. DLA abundances 
apparently can be reduced to a two dimensional problem, as illustrated in Fig. 

4.2. Metal Abundance Ratios 

Metal abundance ratios such as N/Si, Si/Fe, and S/Si are influenced by the effective 
yields of the two metals making up the ratio as well as the relative rates of evolution of the 
stellar groups responsible for their production. Such ratios may also be sensitive to dust 
depletion if one of the elements is more refractory than the other. Problems associated with 
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N/Si and Si/Fe were summarized briefly in the Introduction, where we showed that questions 
surrounding their behavior in DLAs, such as the N/Si scatter and possible bimodality as well 
as the ubiquitous supersolar Si/Fe values, are as yet unresolved. Here we employ the results 
of our models to address some of the questions associated with these abundance ratios. 

4-2.1. N/Si Range and Bimodality 

The large scatter in N/Si described in § 2 is most likely the result of age differences 
among our sample objects. Given the tight link between N/Si and age exhibited by our 
models, we conclude that the scatter in N/Si is evidence of an age range among DLAs, due 
to the delay in N release by LIMS resulting from their slower evolution. Thus, the LNDLAs 
may actually be objects observed at a very young evolutionary age, i.e. a small fraction of a 
Gyr. Initially in these objects the N/Si climbs because of massive star N yields, but as time 
progresses, N produced by LIMS is released during their final evolutionary stages, and N/Si 
rises to greater levels. 

This idea is illustrated in Fig. where we show the time evolution of N following a 
burst of 100 Myr duration. The solid line represents the model with both massive stars and 
LIMS contributing to N buildup, while the short-dashed line shows the track followed by 
a system in which the N production by LIMS has been turned off. The dotted line shows 
the average value for blue compact galaxies on the plateau at -0.28, as inferred by Nava et 
al. (2006) for N/O and converted to N/Si here using the solar O/Si ratio from Asplund et 
al. (2005). The actual range in BCG values extends from -0.39 to -0.12 and is indicated 
with the short vertical line. Recall that the observed range for DLAs of log(N/Si) is -1.16 
to +0.22, shown in the figure by the longer vertical line. For both models, we employed the 
same code as the one described above with SFE=0.01 Gyr -1 . The solid curve turns over at 
roughly the location of the BCG value at a time of 230 Myr and levels off at 400 Myr. HEK 
estimated from their models that most N release has occurred by 250 Myr, although they 
used the LIMS yields of van den Hoek & Groenewegen (1997) which produce more N than 
Marigo's (2001) predictions (see comparison of effective yields of LIMS models by Henry 
2004). In any case it is clear that LIMS contribute significantly to N buildup, according to 
our models and our selection of stellar yields. 

We also conclude from Fig. [7]that BCGs have generally been evolving longer than DLAs, 
since their vertical range is much smaller, and the mean N/Si for the objects on the plateau 
is close to the value at which N/Si of the bursting model above levels off. As already pointed 
out in § 4.1, this conclusion is consistent with other evidence indicating that BCGs are at 
least several Gyr in age. 
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Data for N/Si and Si/H presented in Prochaska et al. (2002) was the first to indicate 
that the distribution of DLAs in the N/Si-Si/H plane may be bimodal. We saw evidence 
of this in Fig. [TJ where a diagonal gap extending from upper right to lower left appears to 
separate the four LNDLAs inside the box from the main body of objects. 

This gap is seen even more clearly in Fig. [HJ where we now plot log(N/Si) versus 
log(N/H). Here the DLA positions are shown with red circles along with the data for blue 
compact galaxies from Nava et al. (2006) indicated with blue crosses. DLA measurements 
which are only limits are shown with green circles and arrows. The four LNDLAs can 
be seen in the lower left, separated from the main group of DLAs by an obvious gap at 
log(N/H) « -7. 

There are two basic hypotheses which would explain such a bimodalitjIE First, it may 
simply be a statistical artifact which occurs when a small number of objects in the N/Si-Si/H 
plane are selected from a large parent population whose members are distributed continuously 
but with non-uniform density within that plane. Or second, the LNDLAs could represent 
a group of objects which evolved with a different set of effective yields than the rest of the 
objects. 

The first hypothesis is related to the explanation originally put forth by Molaro (2003) 
and Molaro et al. (2003). That is, the LNDLAs are objects whose N/Si evolution has 
temporarily been suspended until their LIMS begin releasing N. Once that begins, these 
objects move upward rapidly. Furthermore, with a small sub-population such as ours, a gap 
separating the LNDLAs and the main body of objects might be likely to appear, since the 
probability of finding an object there is small due to rapid evolution through that area. 

We attempted to evaluate this hypothesis by assuming that the parent population of 
DLAs, of which our objects are a subsample, are smoothly distributed in age, i.e. there 
were no times in the past when the formation of systems harboring DLAs was more favored 
than at other times. We calculated five model tracks which differed only in their input SFE. 
These tracks span the entire observed set of objects and are shown in Fig. [HJ The SFE 
values were (upper left to lower right) 0.01, 0.10, 0.50, 1.0, and 4.0 Gyr -1 . Since the track 
corresponding to SFE=0.10 runs through the LNDLAs, we chose it for use in our analysis. 



3 We respectfully acknowledge that some of our colleagues believe that the gap in Fig. [8] is the result of 
small number statistics, that the gap will disappear once more data are available, and therefore the analysis 
which follows is moot. While we certainly agree that this is a possibility we also feel that until that day 
comes, we need to explore the potential explanations for the gap in the event that more data show it to 
be real. If future observations prove that the gap-skeptics are correct, then we'll happily cease our wild 
speculations. 
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Ignoring all objects with upper/lower limits, all DLAs lying close to this track share roughly 
the same SFE but have different ages, where age increases along the track from lower left to 
upper right. We have placed small black circles at intervals of 0.1 Gyr to show how progress 
changes with age. If the DLA population has a uniform distribution of ages, then one would 
expect the same number of objects per interval of spacing of the black circles. 

We have also divided the track into three segments indicated by the four short perpen- 
dicular bars placed along it. The lengths of the segments, beginning in the LNDLA region, 
correspond to 85 Myr, 42 Myr, and 3 Gyr, respectively. There are four objects near the track 
in the first segment, none in the second, and eight in the third. Based on the number of 
objects along the first segment we would expect to find two objects in the gap, where none 
actually occurs. On the other hand, a lack of objects in the gap is very consistent with the 
observation of only eight objects along the 3 Gyr segment. Note that along the 0.10 track 
near the gap there is one object with a N upper limit in the gap and one more like it in the 
direction of greater N/Si. The first point could very well be a LNDLA, while the second one 
could potentially occupy the gap. However, even if this were the case, the results of our test 
would be unchanged. Thus, our test fails to rule out the possibility that the gap is simply a 
statistical artifact. 

We now consider the second hypothesis. In their discussion of the N/Si bimodality, 
Prochaska et al. (2002) used chemical evolution models to demonstrate that this phe- 
nomenon could be explained by a non-standard IMF. They proposed that the LNDLAs 
actually represent a separate population in which star formation is characterized by either 
a highly flattened IMF or an IMF which is truncated below a certain mass threshold. Both 
cases result in a lower proportion of LIMS in the stellar population, and hence the integrated 
effective N yield is reduced. 

In order to check on the feasibility of this second hypothesis, we calculated additional 
individual models for the four LNDLAs in our sample, i.e. Q0151+0448, Q1008+36, Q2206- 
19, and Q2348-14. These models were generally identical to the ones presented in Table [2] 
except that their IMF properties were altered. In the first instance the IMF was truncated 
below 4.80 M Q (or 4.65 M for Q2348-14), while in the second instance the IMF slope was 
assumed to be 0.10 (recall that the IMF slope for the Salpeter relation is 1.35). We then 
proceeded as before in the general modeling, using log(N/Si) and log(Si/H) as our primary 
target values. 

Results of this experiment are presented in Table [31 where the format of Table [2] is 
adopted. Note that the results for the LNDLAs from Table [2J using a normal IMF, have 
been duplicated here for ease of comparison. 
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For the case of a truncated IMF, successful matches to N/Si and Si/H were obtained 
for all four objects. While the IMF for three of the models was truncated below 4.8 M , in 
the case of Q2348-14, whose observed value of N/Si is somewhat higher than the same ratio 
for the other three, the threshold had to be lowered slightly to 4.65 M to provide more N 
to the integrated yield. For the flat IMF, N/Si and Si/H were well-matched for three of the 
LNDLAs but a good model of Q1008+36 proved unattainable. Of the four LNDLAs, this 
object possesses the lowest N/Si value but the highest level of Si/H. Despite numerous trials 
with a wide variety of SFE and age input values we could not find a suitable combination. 

In comparing the SFEs and ages for the three IMF forms represented in Table [3] we see 
that truncated IMF models in general require longer evolutionary ages than either the flat 
or normal IMF examples, since additional time is required to produce the N that would have 
been otherwise synthesized by LIMS. Another trend is that both the truncated and flat IMF 
cases tend to have smaller SFE values. 

We note that Lanfranchi & Matteucci (2003) found that when an altered IMF is used 
to reduce the N/Si to the level of the LNDLAs, Si/Fe becomes unacceptably high. This 
point will be discussed in more detail in the next subsection, where we show that Si/Fe in 
our models does indeed change when the IMF is altered, but not by untenable amounts. 

Therefore, the flattened or truncated IMF hypothesis appears to be capable of explaining 
the LNDLAs. Many theoretical studies of Population III stars clearly predict that the first 
stellar generation either had a greater proportion of massive stars or its IMF was truncated in 
such a way that the relative number of LIMS was small (see the review by Bromm & Larson 
2004). Objects possessing such a stellar population are likely located at high redshifts and 
formed when metallicities were below 10~ 4 Z so i ar (Schneider et al. 2002). Komiya et al. (2007) 
conclude that the changeover from a top-heavy to a normal IMF occurred around [Fe/H] of 
roughly -2.5, i.e. within a dex of the metallicities found for the LNDLAs. Our analysis thus 
suggests that some DLAs may exhibit the chemical imprint expected from Population III 
stars. 



4-2.2. The Supersolar Si/Fe Ratio 

Another important chemical diagnostic of DLAs is the Si/Fe ratio. Since Fe is more 
refractory than Si (Savage & Sembach 1996), the observed gas phase Si/Fe ratio may tend to 
be higher than the intrinsic value. Indeed, DLAs exhibit significantly larger Si/Fe ratios at 
higher Si/H (Figure 2), which is explained by differential depletion (Prochaska & Wolfe 2002). 
For the majority of DLAs in this paper, however, the metallicity is low and the depletion 
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effects may be small. In § 2, we pointed out the tendency for observed Fe abundances to 
be ~ 0.1 dex less than their predicted counterparts in Fig. HI upper right panel. This is 
suggestive of a mild, but non- negligible depletion effect. Note that all of our models predict 
supersolar values for Si/Fe. The average value of Si/Fe predicted by our models is +0.31 
(cr=. 034) for the entire sample, the ratio being roughly constant overall. Fig. [9] is a plot of 
Si/Fe for the models versus the observations. The dashed lines show the solar value (Asplund 
et al. 2005) for each axis. Notice that the observed values are systematically higher than 
the predicted ones and show a much broader range. 

Our models indicate that the supersolar Si/Fe ratios measured in most DLAs stem 
from their relatively unadvanced ages. Roughly half of all Fe is produced in massive stars, 
while the other half comes from Type la supernovae, the precursors of which evolve more 
slowly than massive stars. As a result, Fe release is partially delayed over a timescale of 
between 40-50 Myr for bursting systems or 4-5 Gyr for spiral disks (Matteucci & Recchi 
2001). Therefore, the supersolar Si/Fe values observed in DLAs is expected because of their 
generally young evolutionary ages. On the other hand, the cause of the scatter in Si/Fe 
exhibited by DLAs is more difficult to assess, since both an age range and depletion of Fe 
onto dust are likely to play roles in determining Si/Fe. 

In our discussion of the N/Si ratio above, we pointed out that one objection to the 
IMF modification hypothesis was raised by Lanfranchi & Matteucci (2003), who predicted 
log(Si/Fe) values of roughly +0.6, nearly 0.25 dex above observed levels, for a model with 
an IMF slope of 0.10 while at the same time predicting log(N/Si) levels which were consis- 
tent with those of the LNDLAs. Qualitatively, an increase in Si/Fe is expected, because a 
flattening of the IMF slope reduces the relative contribution of Type la supernovae to the 
effective Fe yield, since these explosive events originate from stars of generally lower mass 
than the stars which produce Si. As a result the Si/Fe ratio should be higher in the case of 
a flattened IMF. 

We investigated this point in some detail, as it may represent a way to eliminate the 
altered IMF hypothesis, thereby clarifying the cause of the N/Si bimodality to some extent. 
We ran models using three different IMF forms, i.e. Salpeter (a=1.35), flattened (a=0.10), 
and truncated (M trunc =4.8 M with a=1.35). For each IMF form we used three different 
values for the SFE, 0.003, 0.05, and 1.0 Gyr -1 . In all cases the calculated value of Si/Fe 
was consistent with the observations, and always remained well below the values predicted 
by Lanfranchi & Matteucci. We speculate that the different results are related either to the 
different yield sets employed or to differences in star formation prescriptions, e.g. bursting 
rather than continuous star formation. We note, however, that regardless of the IMF form 
that we used, we found that values did not differ from observations by more than about 
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0.10 dex and thus within the observational uncertainty of log(Si/Fe). Therefore, we currently 
see no reason to dismiss either the flattened or truncated IMF hypotheses for explaining 
LNDLAs because of a failure to predict observed Si/Fe values. 

In summary, our models employing continuous star formation predict supersolar values 
for Si/Fe that fall within a range which is much narrower than the broad one seen obser- 
vationally. The most likely explanation is Fe depletion onto dust, although bursting star 
formation histories, where burst length and frequency both vary from object to object, may 
also be capable of explaining the observed scatter in Si/Fe. 

4.2.3. S/Si 

The filled circles in Fig. [3] illustrate the tight relation between Si and S observed in 
DLAs. Model results from Table [2] are plotted in the same figure as open circles, where it 
is seen that they all fall very close to the least squares fit (straight line) to the data. Thus, 
observation and theory closely match up in this case. 

The lockstep behavior between Si and S was reported by Prochaska & Wolfe (2002) and 
is consistent with the generally accepted belief that these two elements are forged within the 
same types of stars. In addition, it indicates that Si is not significantly depleted by dust 
formation, at least within the metallicity range of our object sample. The results of Savage 
& Sembach (1996) suggest that while Si can form grains, S generally is not found to be 
depleted in the interstellar medium. If Si depletion were present but unaccounted for in the 
DLA data, then we would expect to find a leftward horizontal offset of the data with respect 
to the models. Clearly, this is not seen. This absence of Si depletion in our data would 
appear to be in conflict with the findings of Dessauges-Zavadsky et al. (2006), who detected 
a positive correlation between [Si/S] and [Zn/Fe], where the second ratio is purported to 
increase directly with dust depletion. The difference may have to do with the fact that we 
consider the total abundances for a set of DLAs whereas Dessauges-Zavadsky et al. (2006) 
examined individual 'clouds' within a DLA. The latter analysis will be more sensitive to the 
depletion effects if there are even modest variations in Si/S and Zn/Fe within a given DLA. 
Furthermore, our sample includes only a few DLAs with [Si/H] > —1 where dust effects 
become more pronounced (Prochaska & Wolfe 2002). 
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5. CONCLUSIONS 

We have compiled published abundances of N, Si, S, and Fe for 45 damped Lyman 
alpha systems in order to study trends present in these cosmologically important objects. 
Our attention was focused particularly on understanding the reasons for: (1) the wide range 
in values of N/Si, their generally low level compared with blue compact galaxies, and the 
apparent bimodality in their distribution; (2) the low metallicity of DLAs as traced by Si; and 
(3) the ubiquitous supersolar Si/Fe values observed. Our primary method of analysis was to 
calculate an individual numerical chemical evolution model for each DLA while attempting to 
keep the number of assumptions about the nature of these objects as well the variable input 
parameters at a minimum. N and Si abundances were our principal targets for modelling, 
and thus we focused our attention on those 30 DLAs for which firm abundance values, as 
opposed to limits, were available. We have reached the following seven conclusions. 

1. The abundances of N, Si, S, and Fe in DLAs are easily explained by a one- zone model 
scheme with infall but no wind, continuous star formation, and only two variable 
parameters, i.e. star formation efficiency and the period of time over which evolution 
occurs (evolutionary age). For now, this 2-dimensional scheme appears to be sufficient, 
and there seems to be no need to increase the level of model complexity. 

2. The ranges in SFE and age for the 30 sample objects modeled are, respectively, 0. OI- 
LS Gyr" 1 and 0.18-2.0 Gyr. The mean of the SFE values is 0.44 Gyr" 1 , while the 
median value for this parameter is 0.20 Gyr -1 . These values agree well with those 
found in previous model studies of dwarf irregulars as well as spiral disk locations at 
large galactocentric distance. 

3. There is a clear anticorrelation between SFE and age for DLAs. Since our models infer 
only the time- averaged SFE, we suggest that star formation in many DLAs may have 
occurred via bursts early on, but the intensity and/or frequency of these episodes has 
subsided with time. In this way older objects will tend to display lower average SFEs. 
Alternatively, the inverse behavior between these two parameters could be due to a 
steady reduction in the interstellar materials necessary for forming stars. 

4. The broad range in N/Si for DLAs compared with blue compact galaxies is explained 
primarily by the fact that we are observing the former at various stages of N enrich- 
ment. This difference between the two object types implies that DLAs are significantly 
younger in evolutionary age than are the blue compact galaxies. 

5. The cause of the apparent bimodality in the N/Si distribution of DLAs is ambiguous, 
as it can be explained either as a statistical result or by the presence of two different 
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sets of effective yields, one associated with a normal IMF and another associated with 
a truncated or flat IMF. Such top-heavy IMFs are qualitatively consistent with current 
theoretical predictions about the mass distributions among Pop III stars, suggesting 
that abundances in the LNDLAs exhibit the imprint of Pop III nucleosynthesis. 

6. The Si/Fe ratio in DLAs exceeds the solar value for all objects in our sample. Thus 
DLAs exhibit a characteristic alpha enrichment also seen in globular clusters and halo 
stars within the Milky Way. However, our models are unsuccessful in explaining the 
large observational scatter in Si/Fe, and we suggest that this may be due to Fe depletion 
onto dust in varying amounts or bursting star formation histories in which both burst 
length and frequency vary among DLAs. 

7. Si and S in DLAs increase together in lockstep, an observational result that is supported 
by our models. Since silicon's abundance is more likely to be depleted through grain 
formation, the observed lockstep behavior strongly suggests that in most DLAs with 
metallicities similar to those in our sample, Si depletion is not important. 

The question of a N/Si bimodality can only be resolved by increasing the size of the 
sample of DLAs with N and Si measurements. More objects will allow us to assess the 
distribution in more detail and hopefully establish whether an apparent bimodality is simply 
the result of a random distribution of a small number of objects or is truly related to 
differences in effective nitrogen yield. Finally, we need to continue to explore the connection 
between DLAs and nearby galaxies. One possible way of searching for such connections 
is to begin with successful chemical evolution models of DLAs and evolve them further in 
time, making predictions both about elemental abundances as well as the integrated spectral 
properties of the more mature galaxy. A comparison of this second generation of models with 
current observations of suspected DLA endpoints such as spiral galaxies and dwarf irregulars 
might reveal an evolutionary connection. 
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A. THE CHEMICAL EVOLUTION MODEL 

Our numerical code is an expanded version of the one employed by HEK. It is a one-zone 
chemical evolution program which follows the buildup of the elements H, He, C, N, 0, Ne, 
Si, S, CI, Ar, and Fe over time. 

We imagine a generic galactic region with a cross-sectional area of 1 pc~ 2 that is initially 
empty but then accretes matter by way of infall of pristine, metal-free gas whose rate is given 
by: 



f (t) = Eto r, 
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where t Q and E to are the current epoch and surface density, the latter is taken to be 5 M Q pc , 
Tscaie = 7 Gyr, and t Q =13.7 Gyifl This form of the infall rate is the one discussed and 
employed by Timmes, Woosley, & Weaver (1995) for their models of the Galactic disk. 

As star formation commences, the total mass M within the region becomes partitioned 
into interstellar gas of mass g and stars of mass s such that 

M = g + s. (A2) 

If ^(t) and e(t) are the rates of star formation and stellar ejection, respectively, then: 

s = m - e(t) (A3) 

and 

g = f(t)-^(t) + e(t). (A4) 
The interstellar mass of element x in the zone is gz x , whose time derivative is: 

gz x + gz x = -z x (t)^(t) + z x f(t) + e x (t), (A5) 

where z x (t) and z x are the mass fractions of x in the gas and in the infalling material, 
respectively, and e x (t) is the stellar ejection rate of x. Solving for z x (t) yields: 

zx(t) = {f(t)[z x - z x (t)] + e x (t) - e(t)z x (t)}g- 1 . (A6) 



4 All models assumed the same surface density of 5 M Q pc 2 , which is consistent with observed values for 
blue compact galaxies, one of the proposed host types for DLAs (L. van Zee, private communication). 
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The second term on the right hand side of eq. A6 accounts for the injection of metals into 
the gas by stars, while the first and third terms account for effects of dilution due to infall 
and ejected stellar gas, respectively. 

The rates of mass ejection e(t) and ejection of element x, e x (t), are: 

-m U p 



e(t) = / [m-w(m)]^(t-r m )0(m)dm (A7) 

•>m Tm 

and 

e x (t) = (1 - C) C up {[m - w(m)]z x (t - r m ) + mp x (m, z t „ Tm )}^(t - r m )0(m) dm 

(A8) 

+C /mw Px ( m ' z t-^) m 0( m ) £ in 24/i 2 ^(t - r m2 )dfi dm B . 

In Eqs. IA7l and lA8l m and z are stellar mass and metallicity, respectively, and m Tm is the turn- 
off mass, i.e. the stellar mass whose main sequence lifetime corresponds to the current age of 
the system. This quantity was determined using results from Schaller et al. (1992). m up is 
the upper stellar mass limit, taken to be 120 M , w(m) is the remnant mass corresponding to 
ZAMS mass m and taken from Yoshii, Tsujimoto, & Nomoto (1996). p x (z) is the stellar yield, 
i.e. the mass fraction of a star of mass m which is converted into element x and ejected, and 
<j)(m) is the initial mass function. Our choice of stellar yields is discussed below. In eq. IA8I 
the first integral gives the contributions to the ejecta of single stars, while ejected masses 
of C, O, Si, S, CI, Ar, and Fe by SNIa through binary star formation are expressed by the 
second integral, where our formulation follows Matteucci & Greggio (1986). Assuming that 
the upper and lower limits for binary mass, Mbm and Mem, are 16 and 3 M , respectively, 
eq. IA8I splits the contributions from this mass range between single and binary stars. The 
relative contributions are controlled by the parameter C which we take to be equal to 0.05 
when 3 < M < 16 M and zero otherwise. The variable /i is the ratio of the secondary star 
mass m2 to the total binary mass, hib, where [i m i n = max(1.5, m Tm , ulb — 8). 

The initial mass function 0(m) is the normalized Salpeter (1955) relation 

1 — a 



m 



(l-a) (1-a) 



m-^ 1+a \ (A9) 



where a=1.35 and Hid OW n = 1 M . Finally, the star formation rate ij)(t) is given by 

m=uM {hT 5 (A10) 

where v is the star formation efficiency in Gyr -1 . The form of eq. A10 is taken from Timmes, 
Woosley, & Weaver (1995). 
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Our calculations assumed a timestep length of one million years, a value which is less 
than the main sequence lifetime of a star with a mass equal to m up , or 120 M . 

At each time point, the increment in z x was calculated by solving eq. A6 and the required 
subordinate equations IA7I and IA8I This increment was then added to the current value and 
the program advanced to the next time step. Finally, the total metallicity at each point was 
taken as the sum of the mass fractions of all elements besides H and He. The unit of time 
is the Gyr, while the mass unit is the solar mass. 

Single staryields for all elements but CI and Ar were taken from Portinari, Chiosi, & 
Bressan (1998)3 for massive stars and Marigo (2001) for low and intermediate mass stars 
(LIMS). An advantage of using these two sets is that they are designed to be complementary. 
The massive star models account for the effects of convective overshooting and quiescent mass 
loss and cover a metallicity range from 0.0004 to 0.05. For metallicities below 0.0004, yields 
were assumed to be equal to their values at 0.0004, i.e. yields were not extrapolated in 
metallicity. However, since N production is exclusively secondary in massive stars, yields 
below 0.0004 were scaled with metallicity. Since Portinari et al. did not publish CI and Ar 
yields, we used the results of Woosley & Weaver (1995). Marigo's computations for LIMS 
account for mass loss, third dredge-up, and hot bottom burning over a metallicity range of 
0.004 to 0.019. Finally, in the case of binary star contributions to evolution of C, O, Si, S, 
CI, Ar, and Fe through SNIa events, we employed the yields of Nomoto et al. (1997), using 
both the W7 and W70 models to account for metallicity effects. 



B. MODEL UNIQUENESS 

In this appendix we consider the issue of model uniqueness within the confines of our 
assumptions, and in particular whether more than one SFE-age combination can adequately 
reproduce the results reported in Table [2j To answer this question, we performed a Monte 
Carlo simulation in which input values of SFE and age were selected randomly from within 
ranges whose extremes matched those of the values inferred for our sample and given in 
Table [21 In this fashion 1000 parameter pairs were constructed. One thousand chemical 
evolution models were then calculated, one for each parameter pair, and the predicted values 
of N/Si and Si/H were plotted along with the observed values in the N/Si-Si/H plane. 

Next, within this plane we selected seven widely spaced DLAs from our sample and 



5 We note that the quantities published by Portinari et al. represent the sum of the amount present in 
the star at birth and that synthesized by the star. This was properly accounted for in our calculations. 
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proceeded to compare our derived values of SFE and age reported in Table [2] with the those 
of nearby models generated by the Monte Carlo simulation. (Since our parameter ranges 
are restricted from the beginning, this is clearly not a foolproof method. However, it should 
at least tell us whether or not there are alternative parameter combinations to the ones in 
Table [2] within those ranges.) 

The results are shown in Table HI where for each of the seven DLAs, we list our derived 
parameters from Table [2] in the first row followed by analogous values for several simulation 
points in subsequent rows. 

The parameter values appear to agree closely for all objects except Q0930+28, where 
our derived age of 1.2 Gyr is more than a factor of two greater than the ages of the nearby 
points. We conclude, however, that in most cases the conditions inferred in the individual 
models above are likely to be unique. 
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Table 1. STATISTICAL SAMPLE: DLA Abundance Summary 



Quasar z aba log N H i log N log Si log S log Fo log(N/Si) log (N/S) log(Si/Fo) Rcf 
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41+° 
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08 
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13 


82 


± 


0. 


02 
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a Calculatcd from log(N/S) and offset by the solar relative abundance, — 0.35dcx. 

Note. — Note that none of the limits reported take into account the uncertainty in the Nhi value. 

References. — 1: Wolfe et al. (1994); 2: Lu ct al. (1996); 3: Prochaska & Wolfe (1996); 4: Prochaska & Wolfe (1997); 5: Lu, Sargent, & Barlow (1999); 6: 
Prochaska & Wolfe (1999); 7: Prochaska & Wolfe (2000); 8: Petitjean, Srianand, & Ledoux (2000); 9: Molaro et al. (2000); 10: Molaro et al. (2001); 11: Prochaska 
ct al. (2001); 12: Ellison et al. (2001); 13: Prochaska et al. (2002); 14: Lopez ct al. (2002); 15: Lcvshakov et al. (2002); 16: Lopez & Ellison (2003); 17: 
Dessauges-Zavadsky ct al. (2004); 18: Prochaska ct al. (2003); 19: Ledoux, Petitjean, & Srianand (2003); 20: O'Mcara ct al. (2006); 21: Dessauges-Zavaksky et al. 
(2006); 22: Prochaska et al. (2007); 23: Lopez et al., in prep; 24: Dessauges-Zavadsky ct al., in prep 



Table 2. CHEMICAL EVOLUTION MODELS 
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log^c/rij 


log(U/±lJ 


l nif /p /u\2 
log^O/rlJ 


10gl_JN/01J 


log^oi/ r c) 


iog(,Wbij 


Q0000-2619 


0.05 


0.50 


-6.68/-6.71 


-6.35/-6.38 


• ■ ■ /-6.72 


-6.66/-6.70 


-4.99/-4.98 


-6.17 


-0.33/-0.33 


+0.31/+0.31 


+0.21 


Q0151+0448 


0.13 


0.20 


-7. 46 /-7. 49 


-6.40/-6.41 


■ ■ • /-6.75 


■ • ■ /-6.76 


■ ■ ■ /-5.00 


-6.21 


-1.06/-1.08 


. . . /+0.35 


+0.20 


CTQ418 


0.25 


0.28 


-6.62/-6.62 


-5.98/-5.97 


-6.51/-6.31 


-6.43/-6.30 


■ ■ ■ /-4.56 


-5.78 


-0.64/-0.65 


+0.45/+0.34 


+0.19 


Q0347-38 


0.18 


0.80 


-S.74/-5.76 


-5.61/-5.62 


<-5.86/-5.96 


-6.13/-5.91 


-4.67/-4.23 


-5.37 


-0.13/-0.14 


+0.52/+0.29 


+0.25 


CTQ247 


0.40 


0.25 


-6.58/-6.55 


-5.81/-5.82 


-6.31/-6.16 


-6.18/-6.16 


. . . /-4.41 


-5.63 


-0.77/-0.73 


+0.37/+0.34 


+0.19 


CTQ247 


0.90 


0.27 


-6.02/-6.01 


-5.50/-5.45 


-5.90/-5.79 


-S.94/-5.78 


■ • • /-4.04 


-5.26 


-0.52/-0.56 


+0.44/+0.33 


+0.19 


Q0450-13 


0.40 


0.22 


-6.77/-6.74 


-5.91/-5.88 


-6.24/-6.22 


-6.38/-6.23 


-5.12/-4.70 


-5.69 


-0.86/-0.86 


+0.47/+0.34 


+0.19 


Q0528-2505 


0.50 


0.26 


-6.37/-6.37 


-5.69/-5.71 


-6.12/-6.05 


-6.01/-6.05 


■ ■ • /-4.30 


-5.52 


-0.68/-0.66 


+0.32/+0.34 


+0.19 


HS0741+4741 


0.10 


0.45 


-6.51/-6.51 


-6.13/-6.14 


-6.47/-6.47 


-6.43/-6.45 


-4.76/-4.73 


-5.94 


-0.38/-0.37 


+0.30/+0.32 


+0.20 


FJ0812+32 


0.12 


0.70 


-6.02/-6.04 


-S.82/-5.85 


■ ■ ■ /-6.19 


-6.il/-6.15 


■ ■ ■ /-4.46 


-5.61 


-0.20/-0.19 


+0.29/+0.29 


+0.24 


FJ0812+32 


0.10 


2.00 


-5.35/-5.46 


-S.37/-5.44 


-5.72/-5.78 


-6.26/-5.68 


-3.86/-4.07 


-5.05 


+0.02/-0.02 


+0.89/+0.23 


+0.39 


Q0841 + 12 


0.50 


0.28 


-6.34/-6.27 


-5.71/-5.67 


-6.18/-6.01 


■ ■ ■ /-6.01 


■ ■ • /-4.27 


-5.48 


-0.63/-0.60 


■ ■ • /+0.33 


+0.19 


Q0930+28 


0.02 


1.20 


-6.55/-6.48 


-6.41/-6.37 


<-6.68/-6.70 


-6.58/-6.64 


-5.01/-4.98 


-6.06 


-0.14/-0.11 


+0.17/+0.27 


+0.31 


Q1008+36 


0.20 


0.19 


-7.40/-7.37 


-6.24/-6.25 


■ • ■ /-6.59 


<-5.59/-6.60 


-4.93/-4.S4 


-6.05 


-1.16/-1.12 


>-0.65/+0.35 


+0.20 


BR1202-07 


0.10 


0.35 


-6.79/-6.76 


-6.25/-6.25 


■ ■ ■ /-6.59 


-6.69/-6.58 


■ ■ ■ /-4.85 


-6.05 


-0.54/-0.51 


+0.44/+0.33 


+0.20 


Q1210+17 


1.30 


0.25 


-S.89/-5.89 


-5.31/-5.33 


-5.64/-5.67 


-5.64/-5.66 


■ ■ ■ /-3.93 


-5.15 


-0.58/-0.56 


+0.33/+0.33 


+0.18 


Q1223+17 


0.23 


0.28 


-6.67/-6.66 


-6.03/-6.00 


■ ■ ■ /-6.34 


-6.34/-6.34 


-6.02/-4.59 


-5.81 


-0.64/-0.66 


+0.31/+0.34 


+0.19 


Q1232+08 


0.40 


0.30 


-6.27/-6.30 


-S.72/-5.74 


-6.07/-6.07 


-6.22/-6.07 


■ ■ ■ /-4.33 


-5.55 


-0.55/-0.56 


+0.50/+0.33 


+0.19 


Q1331 + 17 


0.06 


1.20 


-5.94/-6.00 


-5.86/-5.90 


■ ■ ■ /-6.24 


-6.52/-6.17 


■ ■ ■ /-4.51 


-5.59 


-0.08/-0.10 


+0.66/+0.27 


+0.31 


J1435+5359 


0.20 


0.40 


-6.38/-6.31 


-5.92/-5.90 


-6.27/-6.23 


••■ /-6.22 


■ ■ ■ /-4.49 


-5.70 


-0.46/-0.41 


■ ■ ■ /+0.32 


+0.20 


J1558-0031 3 


0.01 


2.20 


-6.21/-6.33 


-6.43/-6.43 


-6.60/-6.72 


-6.56/-6.61 


-4.81/-5.01 


-6.02 


+0.22/+0.10 


+0.13/+0.24 


+0.41 


J1712+5755 


0.80 


0.18 


-6.71/-6.69 


-S.68/-5.69 


■ ■ ■ /-6.03 


-6.il/-6.04 


■ ■ ■ /-4.27 


-5.50 


-1.03/-1.00 


+0.43/+0.35 


+0.19 


Q1759+75 


1.50 


0.28 


-5.65/-S.67 


-S.23/-5.22 


-5.52/-5.S7 


-5.68/-5.5S 


-4.50/-3.82 


-5.04 


-0.42/-0.45 


+0.45/+0.32 


+0.18 


J2036-0553 


0.05 


0.80 


-6.33/-6.35 


-6.16/-6.17 


■ ■ ■ /-6.50 


-6.52/-6.4S 


■ ■ ■ /-4.77 


-5.92 


-0.17/-0.18 


+0.36/+0.29 


+0.25 


Q2206-19 


0.06 


0.20 


-7.8S/-7.86 


-6.7S/-6.75 


■ • • /-7.09 


-7.il/-7.10 


-S.89/-5.33 


-6.55 


-1.10/-1.11 


+0.36/+0.35 


+0.20 


Q2230+02 


2.00 


0.20 


-5.83/-5.93 


-5.20/-5.26 


-5.56/-5.60 


-5.66/-5.60 


■ ■ ■ /-3.85 


-5.08 


-0.63/-0.67 


+0.46/+0.34 


+0.18 


HE2243-6031 


1.30 


0.28 


-S.79/-5.75 


-5.31/-5.28 


-5.65/-5.62 


-5.75/-5.61 


■ ■ ■ /-3.88 


-5.10 


-0.48/-0.47 


+0.44/+0.33 


+0.18 


Q2342+34 


1.00 


0.22 


-6.21/-6.23 


-5.48/-5.50 


-5.91/-5.84 


-6.08/-5.84 


-5.57/-4.09 


-5.31 


-0.73/-0.73 


+0.60/+0.34 


+0.19 


Q2344+12 


0.10 


0.43 


-6.S7/-6.55 


-6.18/-6.16 


<-6.17/-6.49 


-6.32/-6.48 


-S.33/-4.75 


-5.96 


-0.39/-0.39 


+0.14/+0.32 


+0.20 


Q2348-14 


0.13 


0.23 


-7.21/-7.23 


-6.36/-6.34 


-6.86/-6.68 


-6.74/-6.69 


-5.71/-4.93 


-6.15 


-0.85/-0.89 


+0.38/+0.34 


+0.19 


Average Uncertainty 4 






0.11 


0.089 


0.090 


0.10 


0.27 




0.073 


0.056 




Sun 5 






-4.22 


-4.49 


-4.84 


-4.55 


-3.34 


-3.61 


+0.27 


+0.06 


+0.88 



1 The observed and model values are expressed as observed/model, with typical observational uncertainties of <0.1 dex. 
2 Model results only; No observed values are available. 
3 N yields increased. 

4 Squarc root of the quadrature sum of the logarithmic uncertainties in the observed numerator and denominator values. 
5 Asplund ct al. 2005 



Table 3. CHEMICAL EVOLUTION MODELS: ALTERED IMF 



Quasar SFE (Gyr" 1 ) Age (Gyr) logtN/H) 1 logtSi/H) 1 MS/H) 1 logtFc/H) 1 logtO/H) 1 log(C/H) 2 logCN/Si) 1 logtSi/Fe) 1 log(C/Si) 2 



Truncated IMF 3 



Q0151+0448 


0.006 


3.00 


-7.46/-7.49 


-6.40/-6.46 


• • • /-6.79 ■ ■ ■ /-6.68 


••■ /-5.08 


-6.06 


-1.06/-1.03 ■ ■ ■ /+0.22 


+0.40 


Q1008+36 


0.10 


0.40 


-7.40/-7.32 


-6.24/-6.19 


■•■ /-6.53 <-5.59/-6.51 


-4.93/-4.79 


-5.78 


-1.16/-1.13 -0.65/+0.32 


+0.41 


Q2206-19 


0.003 


3.00 


-7.85/-7.81 


-6.75/-6.76 


■•• /-7.09 -7.il/-6.98 


-5.89/-5.38 


-6.36 


-1.10/-1.05 +0.36/+0.22 


+0.40 


Q2348-14 


0.009 


2.70 


-7.21/-7.13 


-6.36/-6.34 


-6.86/-6.67 -6.74/-6.56 


-5.71/-4.96 


-5.94 


-0.85/-0.79 +0.38/+0.22 


+0.40 


Flat IMF 4 


Q0151+0448 


0.05 


0.20 


-7.46/-7.49 


-6.40/-6.42 


• • • /-6.76 ■ ■ • /-6.88 


••■ /-4.54 


-5.64 


-1.06/-1.07 ■ ■ ■ /+0.46 


+0.78 


Q2206-19 


0.01 


0.40 


-7.85/-7.82 


-6.75/-6.78 


■■■ /-7.13 -7.il/-7.22 


-5.89/-4.91 


-6.01 


-1.10/-1.04 +0.36/+0.44 


+0.77 


Q2348-14 


0.02 


0.51 


-7.21/-7.17 


-6.36/-6.37 


-6.86/-6.71 -6.74/-6.80 


-5.71/-4.51 


-5.60 


-0.85/-0.80 +0.38/+0.43 


+0.77 


Normal IMF 


Q0151+0448 


0.10 


0.20 


-7.46/-7.49 


-6.40/-6.41 


•■■ /-6.75 ■■■ /-6.76 


•■■ /-5.00 


-6.21 


-1.06/-1.08 • ■ ■ /+0.35 


+0.20 


Q1008+36 


0.20 


0.19 


-7.40/-7.37 


-6.24/-6.25 


■•• /-6.59 <-5.59/-6.60 


-4.93/-4.84 


-6.05 


-1.16/-1.12 -0.65/+0.35 


+0.20 


Q2206-19 


0.06 


0.20 


-7.85/-7.86 


-6.75/-6.75 


■•• /-7.09 -7.il/-7.10 


-5.89/-5.33 


-6.55 


-1.10/-1.11 +0.36/+0.35 


+0.20 


Q2348-14 


0.13 


0.23 


-7.21/-7.23 


-6.36/-6.34 


-6.86/-6.68 -6.74/-6.69 


-5.71/-4.93 


-6.15 


-0.85/-0.89 +0.38/+0.34 


+0.19 



1 The observed and model values are expressed as observed/model. 
2 Model results only; No observed values are available. 

3 IMF is truncated below 4.80 M© for all but Q2348-14, where 4.65 M is used. 
4 The IMF slope is a=0.10 
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Table 4. MODEL UNIQUENESS TEST RESULTS 



Data Source 1 SFE (Gyr -1 ) Age (Gyr) 



Q2206 



Model 


0.060 


0.20 


SI 


0.046 


0.23 


S2 


0.039 


0.28 


S3 


0.065 


0.21 


S4 


0.039 


0.32 


S5 


0.070 


0.20 


Q0151 


Model 


0.13 


0.20 


SI 


0.14 


0.20 


S2 


0.10 


0.25 


S3 


0.094 


0.26 


S4 


0.13 


0.19 


Q0930 


Model 


0.02 


1.2 


SI 


0.050 


0.49 


S2 


0.058 


0.42 


S3 


0.05 


0.49 


Q1759 


Model 


1.50 


0.28 


SI 


1.36 


0.30 


S2 


1.33 


0.31 


S3 


1.32 


0.31 


S4 


1.28 


0.32 


S5 


1.34 


0.30 


Q0450 


Model 


0.40 


0.22 


SI 


0.39 


0.21 


S2 


0.38 


0.21 


S3 


0.35 


0.23 


S4 


0.44 


0.20 


S5 


0.34 


0.23 


CTQ247 (z aba =2.595) 


Model 


0.90 


0.27 


SI 


0.72 


0.30 


S2 


0.78 


0.28 


S3 


0.81 


0.27 


S4 


0.78 


0.27 


S5 


0.66 


0.31 


Q1210 


Model 


1.30 


0.25 


SI 


1.39 


0.25 


S2 


1.39 


0.25 


S3 


1.33 


0.26 


S4 


1.28 


0.27 


S5 


1.40 


0.25 



1 Values of SFE and age pertaining to the 
model in Tablc[2]arc listed in the first line 
for each object. Each subsequent line pro- 
vides values for several simulation points. 
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Fig. 1. — log(N/Si) versus log(Si/H) showing the positions and uncertainties for DLAs in 
our sample with firm abundance measurements (filled circles), as well as those objects with 
upper (down triangles) or lower (up triangles) limits. Also included are the measurements 
for a sample of blue compact galaxies from Nava et al. (2006) in which their O abundances 
have been transformed into Si abundances using the O/Si solar abundance ratio of Asplund 
et al. (2005). The double-arrowed vertical line indicates the extent of the N/O plateau for 
BCGs as defined by Nava et al. The box in the lower left region surrounds the low nitrogen 
DLAs. 
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Fig. 2. — log(Si/Fe) versus log(Si/H) showing the positions and uncertainties for DLAs in our 
sample. Globular cluster data from Pritzl et al. (2005) are indicated with crosses. The solar 
value for log(Si/H) is -4.49, outside the range of the plot. The horizontal dashed line shows 
the solar abundance ratio for log(Si/Fe), where the solar values were taken from Asplund et 
al. (2005). 
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Fig. 3. — log(S/H) versus log(Si/H) showing observed values and uncertainties for our sample 
objects as filled circles. Model-predicted values for the same objects are shown as open circles. 
The bold line shows the results of a least squares fit to the observations, with values (and 
uncertainties) for the y- intercept (A), slope (B), and correlation coefficient (r) given in the 
upper left corner. 
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Fig. 4. — Residuals between the model predictions and the observed values versus observed 
values for log(N/H), log(Si/H), log(Fe/H), and log(S/H). The absolute value of the average 
residual is indicated in each panel. The horizontal dashed lines show the zero residual levels. 
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Fig. 5. — Model-determined values for the star formation efficiency in Gyr 1 and evolutionary 
age in Gyr for our sample DLAs. A box has been drawn around the low nitrogen DLAs. 
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Fig. 6. — Upper panels: Observed values of log(N/Si) versus SFE, age, and the product of 
these two parameters, respectively left to right. Lower panels: Same as the upper panels but 
for log(Si/H). 
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Fig. 7. — log(N/Si) versus the time in Gyr following a starburst lasting 100 Myr. The solid 
line relates to a model in which both massive stars and LIMS contribute to N production, 
while the short dashed line shows the result for a similar model except that N production 
by LIMS is turned off. The long horizontal dotted line indicates the average log(N/Si) for 
plateau blue compact galaxies, using log(N/0) from Nava et al. (2006) and converting to 
log(N/Si) using the solar values in Asplund et al. (2005). The shorter of the two vertical 
lines on the right shows the BCG range, while the longer one shows the range for DLAs. 
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Fig. 8. — Observed values of log(N/Si) versus log(N/H) for our sample of DLAs (red circles) 
and the blue compact galaxies (blue crosses) from Nava et al. (2006). DLA values where 
only limits are available are shown with green dots with arrows. The five model tracks 
show predicted evolution through the plane for different star formation efficiencies, where 
the values of the latter in Gyr -1 , upper left to lower right, are 0.01, 0.10, 0.50, 1.0. and 4.0. 
Time steps in increments of 0.1 Gyr are shown with black circles on the track corresponding 
to SFE=0.10. This same track is divided into three segments by the short bars perpendicular 
to the track (see text). 
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Fig. 9. — log(Si/Fe) values from models (vertical axis) versus observed values (horizontal 
axis). The dashed lines show solar values in each direction, from Asplund et al. (2005). 



